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Abstract. We study a dephasing channel with memory, described by a Hamiltonian 
model in which the system-environment interaction is described by a stochastic process. 
We propose a useful way to describe the channel uses correlations. Moreover, we give 
a general expression for the coherences decay factors as a function of the number of 
channel uses and of the stochastic process power spectrum. We also study the impact 
of memory on the three qubit code, showing that correlations among channel uses 
affect very little the code performance. 



1. Introduction 

State transfer between different units of a quantum computer or entanglement 
distribution between two parties require quantum communication channels [H [2]. 
They are quantum systems transferring quantum information: the proper quantity to 
characterize the channel performance is the quantum capacity, defined as the maximum 
number of qubits that can be reliably transmitted per channel use [3]. 

Quantum channels are often thought as memoryless, implying that the effect 
of the channel on each information carrier is always described by the same map 
S. In other terms there is no memory in the interaction between carriers and the 
environmental degrees of freedom physically describing the channel. In this case the 
quantum operation for channel uses is given by Sn = . However, in several 
physically relevant situations this is not a realistic assumption. Memory effects appear 
when the characteristic time scales for the environment dynamics are comparable or 
longer than the time between consecutive channel uses. For instance, solid state 
implementations, which are the most promising for their scalability and integrability, 
suffer from low frequency noise [1]. In optical fibers, memory effects may appear 
due to slow birefringence fluctuations [5]. This introduces correlation among uses, 
then £m 7^ this kind of channels being referred in the literature as memory 

channels O [71 [8] . 



Semiclassical model for a memory dephasing channel 



2 



A very interesting question, raised for the first time in Ref. [H], is wlietlier 
memory can enhance the transmission capacity of a quantum channel Recently we 
have considered a channel subject to dephasing noise described by a Markov chain, 
showing that the quantum capacity increases with respect to memoryless limit [10] (see 
also pn fT2]). Furthermore based on theoretical arguments and numerical simulations, 
we have conjectured that the enhancement of the quantum capacity also takes place for 
a dephasing quantum environment modelled by a bosonic bath |10j . 

This issue is also relevant for the performance of Quantum Error-Correcting Codes 
(QECCs). Since quantum capacity is the maximum rate of reliable quantum information 
transmission, it puts an upper bound to the asymptotic rate achievable by any QECC. 
On the other hand, realistic QECCs necessarily work on a finite number of channel 
uses. Moreover, present day experimental implementations [13l [14j are bases on very 
few channel uses. Previous studies have investigated the impact of correlations on the 
performance of QECCs [151 [16]. Depending on the chosen model, correlations may have 
positive or negative impact on QECCs. In a previous paper [17] we have shown, for a 
Markovian dephasing channel, that also low values of memory, for which the quantum 
capacity does not change appreciably, can have a detrimental impact on the three-qubit 
code performance. 

In this paper we describe a dephasing channel by a Hamiltonian where the system 
environment interaction is modelled by a stochastic process. Then we discuss the three- 
qubit code error performance in presence of channel correlations. 

2. Channel Model 

We suppose that information is carried by qubits that transit across a communication 
channel, modelled as an environment determining pure dephasing of the qubits. The 
environment acts as a stochastic drive ^{t) on the system and the Hamiltonian describing 
the transmission of qubits through the channel reads 



The k—th qubit is coupled to the environment via its Pauli operator al , with coupling 
strength A. The functions fk(t) = u{t — t^) — u{t — — Tp), where u(t) is the unit step 
function [18], switch the coupling on and off. Here Tp is the time each carrier takes 
to cross the channel; r = t^+i — tk is the time interval that separates two consecutive 
qubits entering the channel. Only when the A;-th qubit is inside the channel the function 
fk = I- We assume ^(t) is a stationary and Gaussian stochastic process [19] with zero 
average value, characterized by its autocorrelation function C(r). 

To deal with this problem, we first consider the time evolution of the system for a 
given realization ^{t) of the stochastic process, and then we perform an average over all 




(1) 



k=l 
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possible realizations. The A^-qubit time evolution operator for a given realization is 

N 

U^{t) = (g)exp(-iaf 0;.), (2) 

k=l 

where (pk is the phase acquired by the k-th qubit coherences after the qubit crossed the 
channel: 

0, = ^y^ at')dt'. (3) 

Time evolution is conveniently described in the factorized basis {\j) = 
jjv), ji, JTV = 0,1}, where {\jk)} are eigenvectors of Let = 

Ylj be the initial state of the A^-qubit system; the final state after all 

TV qubits crossed the channel is given by 

pf = U^it)pQulit) = (21 (4) 

j,l ^ k=l ' 

where = Ik — jk = ~ By averaging over the stochastic process we 

finally obtain 

N 

p""' = (pf) = ya,i ( exp ( 2i Vsk^k ] ) (5) 



^' = (pf ) = Yl \ ( 5Z ^^^^ ] 
j,i ^ ^ k=i ^ 



which is a quantum operation for the N-qubits system: p®' = E^{p'^\ It is possible to 
show that the quantity XlfcLi ■^a.-^a,- is itself a Gaussian variabl^, therefore 

/ exp ( 2i ^ Sfc^fc ) ) = exp f - 2 ^ SkSk'{(t)k(i>k') J • (6) 
\ ^ fc=i ^ ' ^ k,k'=i ^ 

We call this quantity the (j, /)-coherence decay factor, since it is just the damping 

experienced by the (j, /) system coherence: 

Next, by using the stationarity of ^(t), we calculate 

{<Pk<Pk') = {-r dh^ti) / dt2at2: 

\ 4 Jtk Jty 
— dti / dt2 C(ti-t2). (8) 

Since the autocorrelation function can be expressed in terms of the power spectral 
density S{uj) = J dre^'^^C(r) of ^(t) we obtain 

, V .9 f°° ^, , 1 — COsfcUTp) r /, ,/x 1 

{(pk(pk') = / — cos cu(A; - k')T], 9 

Jo 27r u;^ 

if Given two Gaussian variables x and y with arbitrary variances and some degree of correlation, the 
two variables z± = x ± y are again Gaussian, as one can find by deriving the (z+ , z_ ) mixed density 
function from the one of {x,y) [19] . Each phase 0^ is a time integral of a Gaussian stochastic process, 
so it can be view as the limit of a sum of Gaussian variables, so in its turn it is Gaussian. 
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and the final form of the coherences decay factor follows: 



oo 



Dji = exp - AM — S{uj) 





(4 - jfc)(/fc' - jfc') cos[uj{k - k')T] j (10) 

k,k'=l ^ 

This result is identical to the coherences decay due to a dephasing channel modelled by 
a set of quantum harmonic oscillators [10] . 

The expression ffTOl) can be put in a nice and useful form. To this end we define 
/ifcfc/ as the correlation coefficient [I9] between the phases 0^ and 0fc/: 

where we set r/^ = {(pl); in fact, thanks to stationarity of the process, the quantity 
does not depend on k (see equation (Q). The coherence decay factor can be 
rewritten as 

N 

Dji = exp - 2r]^ ^ SkSy IJ^kk') ■ (12) 

Now we observe that g = exp(— 2?]^) is just the damping experienced by single qubit 
coherences for one channel use (A^ = !)• We finally write 

where we have defined ^ik-w = fJ^kk'- In fact the stationarity of ,^(if:) imphes that fikk' 
depends only on \k — k'\. The quantity fik-k' is a measure of the degree of the correlation 
between the channel uses k and k'. 



3. Three-Qubit Code performance 

As a measure of the quantum information transmission reliability we use the 
entanglement fidelity [2U]. To define this quantity we look at the system Q as a part 
of a larger quantum system TZQ, initially in a pure entangled state The initial 

density operator of the system Q is then obtained from that of TZQ by a partial trace 
over the reference system TZ: = Ttti[\iIj'^^) {ip'^^l]. The system Q is sent through the 
channel, while TZ remains ideally isolated from any environment, being p^®' the final 
state of TZQ after the transmission. Entanglement fidelity is just the fidelity between 
the initial and the final state of TZQ: 

First we consider a single use of the channel described by Hamiltonian ([1]). We suppose to 
feed the channel with a quantum source [21] described by the density operator = jl. 
The entanglement fidelity is |22j : 
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where I''^ is the identity operator and = Si. This case is relevant in quantum 
information field as it takes place when two communication parties try to share a Bell 
state: the party that initially possesses the pair sends one half of it through the quantum 
channel S. is the fidelity between the actually shared pair and the original one; it 
means that a Bell measurement on p^®' able to distinguish the ideally shared state from 
the other states of the Bell basis fails with error probability = 1 — F^. From ([T5j) 



it follows that the error probability for a single channel use is in what follows we 
identify this quantity by e. 

Now we suppose to use the Three-Qubit Code (TQC) [Tl[2] to send p®. The system's 
state is encoded by using two ancillary systems A and B. The system and the ancillary 
qubits are encoded by means of a set of quantum operations that we resume as C^-^^ 
(stages a, b, c in figure [1]) and then transmitted in iV = 3 uses of channel ([T]). Then the 
receiver performs the decoding (stages e, f, g, h in figured]) on the system QAB. 

After tracing out AB, he obtains the final, generally mixed state of system TZQ: 

where IV^^^^^) = I^A^^) ® |00-^^). Entanglement fidelity fF^^'^ = {^^'^\Ptqc\^^^) 
just gives the probability that the code is successful. The merit of this code is that it 
drastically reduces - in absence of use correlations, i.e. for £Q-^^ = Sf'^ - the transmission 
error probability from e to Pf 
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Figure 1. Scheme of a three qubit code [13]. This quantum error correcting code 
was initially designed for a bit flip channel, for this reason each channel use (stage 
d) is embedded between two Hadamard gates [HE] (stages c and c). The coding is 
performed by means of CNOT gates (stages a, b, f and g), the decoding also requiring 
a Toffoh gate [lll2] (stage h). 



Now we investigate the effects of channel correlations on the performance of a TQC. 
After some involved calculations it comes out that 

1 3 1 



2 4^ \<o 



3\ n'^fJ-QA-^fJ.QB-'^IJ-AB _| 



g-'ifJ-QA+'ifJ-QB-^HAB _j_ 
-I'ifJ.QA+'^lJ-QB+'^lJ-AB] 



By observing that /iQ_4 = /i^g = pi and pge = p2 we can rewrite equation f[T7|) as: 

1 3 1 



2 4 16^ 
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An interesting result can be obtained by considering the case of a small error probability 
e ^ 1. In this regime we can take the series expansion of f|T8|) near e = 0: 

piTQCm) ^ 1 _ (3 + ^ 2i2l)e'. (19) 

This expression tells us that even though memory lowers the fidelity, this worsening is 
always slight and absolutely negligible when yUi,/i2 ^ 1- Moreover, it highlights that 
channel correlations - inside the Hamiltonian model ([T]) - permit the TQC to maintain 
its error probability p^'^^'^''^^ = i _ p('^Q'-^'"^) Qf ^j-^g order of e^. However, one has to 
take care that in the case perfect memory the code error triplicates from 3e^ to 9e^. 
These results are very similar to the ones by Clemens et al. [16]. However, we discuss 
time rather than space correlations and average with respect to stochastic processes. 

Rather than choosing a particular autocorrelation function C(r) for ^(t) and then 
trying to carry out a specific relation between it and the TQC error probability, we 
make some general considerations about the impact of correlations on the code error 
probability. For the sake of simplicity we assume /ii and /i2 having positive values 
(we do not consider anti-correlation cases). While the range of /ii is [0,1], we can 
argue that /i2 < /ii since one expects that the phase correlation does not increase 
when increasing the channel uses distance; furthermore it can be proved that it must 



be /i2 > /i2 = — ly. Studying the first derivatives of p^'^'^'-'"'"^^ as a function 
of Hi and /i2 it turns out that p^'^^^'"^^ is monotonical with respect to /ii (error 
grows with /ii), but not with respect to fi2- It indeed can displays a minimum at 
/i2„pj = —0.25\ogg[{g^^^ + g~'^'^^)/2], but its presence is substantially irrelevant, and 
one can reasonably say that the code error probability is also increasing with respect 
to /i2. Thus to characterize the p^^^^'"^^ behaviour, we plot it as a function /xi, using 
/i2 G {max(0, /i2), /ii} as parameter. As it is showed in figure[2l in which we set e = 10~^, 
the TQC error probability weakly depends on /ii, and the ^2 allowable values affect very 
little it. In the same figure we also plot the error probability for a two-qubit code [22] 
encoding a qubit into the subspace spanned by {|01), |10)}: the performance of this last 
code is always worse than the TQC ones, unless in the case /xi — 1, for which the coding 
subspace becomes decoherence-free. 

The TQC exhibits the same kind of behaviour showed in figure [2] as e changes. In 
figure [3] we plot p^^^'^c.m) ^ function of e for /i2 = yUi = 1, the case in which the code 
shows the worst performance; we do not plot the cases of low correlations (/ii < 0.1) 
since the correspondent curves are practically indistinguishable from the memoryless 
ones. We also compare p^^^'^'^^ with the error probability of the two-qubit code [22j : 
to produce good results this last one requires very high degrees of correlation between 
successive channel uses. 

In conclusion we find that a Hamiltonian formulation of a memory dephasing 
channel shows that the three qubit code is robust against channel correlations. 

§ One can see it by considering the average of [02 + a(0i + '/'s)]^ where a is a real variable: it is a 
quadratic form in a and by imposing that it must always be positive we obtain the desired condition 
on fi2 [El- 
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Figure 2. Plot of code error probability as a function of fii . The dotted and the dashed 
grey lines represent respectively the single channel use error probability (e = 10"'^) and 



TQC error probability for the memoryless channel (Pc 



The error probabilities 



for the three qubit code in presence of correlations (p^QC".™) are represented by black 
curves: solid curve refers to = rnax(0, /22) and the dotted one to H2 = fJ-i- There 
is also displayed (solid gray curve) the error probability for a simple two-qubit code 
encoding a qubit into the subspace spanned by {|01), |10)}. 



Significantly different results emerge if we describe channel correlations inside a 
Markovian model [I7j: In this latter case memory restores the e— dependence in the 
code error probability, thus drastically reducing the code performance. 



ie+00 



ie-02 



le-04 



Oh 



le-06 



le-08 
le 



04 



le-03 



le-02 

e 



ie-01 



Figure 3. Plot of code error probability as a function of the single channel use 
error probability e. The dashed grey line represents the TQC error probability in 
the memoryless case (Pe'^^"^. For the error probabilities of the three qubit code 
in presence of correlations (pj^-^'^'^'™) ) -^^g p^Q^ tj^e worst case (/X2 = /ii = 1) by a 
dotted black curve. There is also displayed (solid gray curve) the error probability for 
a simple two-qubit code encoding a qubit into the subspace spanned by {|01), |10)} for 
fii = 0.99 (triangles down). 
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